Evolution of the parton dihadron fragmentation functions 
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Quark and gluon parton dihadron fragmentation functions and their evolution are studied in 
the process of e'^e" annihilation. We provide definitions of such dihadron fragmentation functions 
in terms of parton matrix elements and derive the momentum sum rules and their connection to 
single hadron fragmentation functions. We parameterize results from the Lund Monte Carlo model 
JETSET as the initial conditions for the parton dihadron fragmentation functions at the scale 
Qo = 2 GeV^. The evolution equations for the quark and gluon fragmentation functions are solved 
numerically and the results at different higher scales agree well with JETSET results. The 
importance of the input from the single fragmentation functions is pointed out. 
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I. INTRODUCTION 

In the analysis of particle production from jets in 
e"'"e~annihilation or in p + p{p) collisions, a primary 
observable is the single particle inclusive cross sec- 
tion. Within perturbative QCD (pQCD) and at leading- 
twist, the single particle inclusive cross section in 
e"'"e~ annihilation can be proved to factorize into a cal- 
culable perturbative hard partonic part and a nonper- 
turbative single inclusive fragmentation function In 
elementary particle collisions it is also possible to mea- 
sure multiple particle production cross sections. This has 
led to the construction and study of multi-particle ob- 
servables such as event shapes 0il3- These offer further 
insight into the substructure of jets produced in high en- 
ergy particle collisions. In the analysis of jets produced 
in high energy heavy-ion collisions or in semi-inclusive 
deeply inelastic scattering (DIS) off large nuclei, however, 
the latter analysis is quite infeasible. Along with the 
single inclusive cross section another measurable quan- 
tity which may offer insight into the modification of jet 
properties in a medium are two high momentum particle 
correlations. 

Dihadron correlations have indeed been measured re- 
cently in a variety of experiments. In the study of jet 
suppression in heavy-ion collisions, correlations between 
two high pt hadrons in azimuthal angle arc used to study 
the medium modification of jet structure in heavy-ion 
collisions at the Relativistic Heavy-ion Collider (RHIC) 
0,13 . While the back-to-back correlations are suppressed 
in central Au + Au collisions, indicating large parton en- 
ergy loss in the dense medium, the same-side correlations 
remain approximately the same as in p + p and d + Au 
collisions. This is to be contrasted with the large suppres- 
sion observed in single inclusive spectra . Two hadron 
correlations within the same jet are also measured in DIS 
off various nuclei at the HERMES experiment 0. In 
such experiments the dihadron correlations surprisingly 
display minimal variation with the choice of nuclear tar- 
get, even though the single inclusive production of lead- 
ing hadrons are significantly suppressed with increasing 



atomic number of the nuclear target Given the ex- 
perimental kinematics of both experiments, this may be 
considered as an indication of parton hadronization out- 
side the medium. However, since the same-side correla- 
tion corresponds to the two-hadron distribution within a 
single jet, the observed phenomenon is highly nontrivial. 

To study systematically such phenomena, an extension 
of the single inclusive fragmentation formalism of QCD 
is required, to include correlations between pairs of par- 
ticles produced in the same jet. Toward this end, we had 
proposed dihadron fragmentation functions in a recent 
paper 0. In that effort, the double differential cross 
section for the same-side production of two hadrons in 
the e"'"e~ annihilation was factorized into the well known 
hard cross section for e+ + ^ 7* ^ qq and the quark 
dihadron fragmentation functions. The essential purpose 
of this fragmentation function is to measure the distribu- 
tion of hadron pairs produced in the fragmentation of a 
hard quark. 

As shown in the case of single hadron fragmentation 
functions [Tol |. the medium modification of the fragmen- 
tation function due to multiple scattering and induced 
gluon radiation closely resembles that of radiative cor- 
rections due to gluon bremsstrahlung in vacuum. There- 
fore, it is important to understand first the QCD evo- 
lution of the dihadron correlations in quark and gluon 
jets in vacuum. In the last study 0, we had concen- 
trated on the definition of the dihadron fragmentation 
functions and derived the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) evolution equations [Tlllllll^ 
for the non-singlet quark dihadron fragmentation func- 
tion. In this paper, we will consider singlet quark and 
gluon dihadron fragmentation functions which will be di- 
rectly relevant to the study of medium modification of 
jets in high-energy heavy-ion and DIS collisions. In ad- 
dition, we will explore momentum sum rules for the di- 
hadron fragmentation functions and their relationship to 
single hadron fragmentation functions. 

In the factorized form for the inclusive hadron pro- 
duction cross section, the short distance parton cross 
section can be computed order by order as a series 
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in Q!<,(Q^) for reactions with momentum transfer much 
higher above Aqcd- The long distance objects or the 
inclusive n-hadron fragmentation functions contain the 
non-perturbative information of parton hadronization 
P, i3| . These fragmentation functions can be defined in 
an operator formalism 0, ^| and hence are valid be- 
yond the perturbative theory. They however cannot be 
calculated perturbativcly and have to be instead inferred 
from experiments. Since the definitions of these functions 
are universal or process independent, once measured in 
one process, e.g. e+e~ annihilation, they can be applied 
to another, e.g. deep inelastic scattering or p + p col- 
lisions, and therein lies the predictive power of pQCD. 
Yet another predictive power of pQCD rests in the fact 
that once these fragmentation functions are measured or 
given at one energy scale, they can be predicted for all 
other energy scales via the DGLAP evolution equations. 



While hadronization of the out going quark and gluon 
jets is a non-perturbative phenomena, inclusive hadron 
production cross sections in e+e~ collisions have turned 
out to be one of the many successful predictions of per- 
turbative QCD [l3,|ll|l3|. To the best of our knowledge, 
however, measurements of dihadron fragmentation func- 
tions have not been performed in e^e~experiments. In 
the numerical study of the nonsinglet quark dihadron 
fragmentation functions in Ref. |9j, we used a simple 
ansatz for the initial condition as Dpfsizi, Z2) = D{zi) x 
D{z2), which at best was just a guess and, as we will 
show later, differ significantly from the inherent hadron 
correlations in a single jet. To facilitate a more accurate 
numerical study of the singlet quark and gluon dihadron 
fragmentation functions in this paper, we will parameter- 
ize the results of the Lund Monte Carlo model JETSET 
[Tgj which has successfully described many aspects of jet 
fragmentation in both e"'"e~and DIS experiments. 



The remaining sections of this paper are organized as 
follows. In Sec. II we review the definition of the double 
hadron fragmentation function. We outline how such a 
function may be isolated in the expression for a double 
differential inclusive cross section. In Sec. Ill we derive 
and discuss various sum rules that are obeyed by the di- 
hadron fragmentation functions and their connection to 
single hadron fragmentation functions. In Sec. IV we 
outline the derivation of the DGLAP [H [Hill evolu- 
tion equations for the quark and gluon dihadron fragmen- 
tation functions. The initial conditions for the dihadron 
fragmentation functions are extracted from JETSET at 
a scale Qo in Sec. V. They are then evolved numerically 
via the DGLAP evolution equations to different scales Q 
and compared to the result form JETSET. Finally in Sec. 
VI we discuss the results of our calculation and present 
our conclusions. 



II. DIHADRON FRAGMENTATION 
FUNCTIONS 

In this section, wc will review the definition and prop- 
erties of dihadron fragmentation functions in the the 
semi-inclusive process e+ -I- e~ 7* ^ hi + h2 + X 
via single jet fragmentation where two hadrons are iden- 
tified. In most cases we will be concerned with back- 
to-back quark and antiquark jets. Study of the gluon 
fragmentation function, however, will involve three jet 
events. In both cases, our focus will always be on two 
hadrons produced in the same jet. 

The cross section for the process of e"'"e~annihilation 
into hadrons may be expressed as 

Shad 

X C^.iOlJ^'mShad) {Shad I J" (0) 1 0) 

~ 2s(74 4 ' ^ > 

where = tpqj'^ipq is the hadronic electromagnetic 
current; C^i, the leptonic tensor and W^"^ the hadronic 
tensor. The four momentum of the virtual photon is 
q = ki + k2 = ((5,0,0,0) and the Mandelstam variable 
s ~ = Q'^. In the remaining, the sum over Shad will 
include both the sum over the complete set of states and 
the phase space integration n/g5^^^d'^p//2i5/(27r)'^ and 
Ps = EfP/- 



Pi 




FIG. 1: The leading order Feynman diagram contributing to 
the double inclusive fragmentation function. 

The definition and factorization of dihadron fragmen- 
tation functions involve identifying two hadrons with 
nearly parallel momenta pi and p2 among hadronic states 
along the direction of one of the partons and replacing 
the remaining sum over hadronic states with a sum over 
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the rest of all partonic states (see Fig. ^ . This is fol- 
lowed by an extraction of the leading twist component 
(see Ref. for details). Differentiating the total inclu- 
sive cross section with respect to the forward momentum 
fractions of the two hadrons i.e., 2:1,^2, the leading or- 
der (LO) factorizcd form of the double differential cross 
section may be expressed as 



dzidz2 



D. 



hih'2 



{Z1,Z2) 



(2) 



where ctq"^ is the total LO hard cross section for an 
e"'"e~pair to annihilate into hadrons. The two functions 
0^1^2(2^,22) and £'gi''2(zi, Z2) represent the dihadron 
fragmentation functions for a quark and an antiquark. 
The quark dihadron fragmentation function in light-cone 
gauge (A+ = gauge) is obtained as, 




D 



hi 



(21,22) 



4,ziZ2 J 4(27r)3 J (27r)4 
X Tq(j);pi,p2)6 (^Zh - ^ 



(3) 



The forward momentum fractions of the identified 
hadrons zi, Z2 are zi = p^ /p^ and Z2 = P2 jp^ ■ The mo- 
mentum ph represents the sum of the hadronic momenta 
i.e., ph = Pi + P2- Henceforth, the direction identified 
by Ph will be considered the same as the direction of the 
jet. The transverse spread of the two hadrons around 
the jet direction is indicated by the relative transverse 
momentum q± = pi^ — P2±- As is indicated by the final 
(5-function in the above equation, the momentum fraction 
Zh = 2:1-1-2:2. The matrix element of the operator product 
Tq(p; 2i, Z2) is given as (in the light-cone gauge) 



Tq{p;pi,P2) = Tr 



x{0\ijq{x)\pi,p2, S - 2){p,,p2, S - 2\i;q{0)\0) 



(4) 



The above expression for the dihadron fragmentation 
functions includes the integration over the transverse mo- 
menta q^. As a result the angular correlation between 
the detected hadrons is also integrated over in such a 
definition. In order to observe such correlation a dif- 
ferentiation of the dihadron fragmentation function with 
respect to the transverse angle must be carried out. In 
this paper, we will continue to focus on the integrated 
fragmentation function. 

The above definition for the dihadron fragmentation 
function also admits a simple diagrammatic interpreta- 
tion in terms of the cut- vertex method of Mueller [l^ . 
The diagrammatic rule in terms of the quark cut-vertex is 
outlined in Fig. [3 The product of these rules along with 
the factor z^/ (4ziZ2) is the integrated dihadron fragmen- 
tation function. 



FIG. 2: The cut-vertex representation of the quark dihadron 
fragmentation function. 



The gluon dihadron fragmentation function can be 
similarly constructed by identifying two hadrons mov- 
ing with almost parallel momenta in the direction of the 
outgoing gluon in a three jet event in e"'"e~ annihilation 
processes. Factorizing the hard cross section from the 
soft matrix element one obtains the gluon dihadron frag- 
mentation function at leading twist as (see Ref. |^ for 
details). 



Dg{Zl,Z2) 



dH 

2ziZ2 J 8(27r)2 J (27r)4 
X S{zh-^)Tg{l;pi,p2). (5) 



In the above equation, the meanings of various momenta 
and forward momentum fractions arc the same as for 
the quark dihadron fragmentation function. The gluon 
overlap matrix element Tg(/;pi,p2) is defined as 

Tg{l;pi,P2) = j d''xe' ''Y.^0\Al{x)\p,,p2,S-2) 



s-2 

X (pi,p2,5-2|^^(0)|0) 



16 ' 



(6) 



where d^^{l) is the gluon's polarization tensor in the 
light-cone gauge and sum over the color indices of the 
gluon field is implied. 

The gluon dihadron fragmentation function also ad- 
mits a simple diagrammatic interpretation in terms of 
the cut-vertex method. The diagrammatic rule in terms 
of the gluon cut-vertex is outlined in Fig. |21 The prod- 
uct of these rules along with the factor z\l(2zxZ2) is the 
integrated dihadron fragmentation function. 
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FIG. 3: The cut- vertex representation of the gluon dihadron 
fragmentation function. 



III. SUM RULES 

One of the many interesting properties obeyed by the 
single inclusive fragmentation functions are the various 
sum rules. Primary among them are the momentum sum 
rules 



(7) 



Here the sin gle inclusive fragmentation functions are de- 
fined as lilMEl, 



(2^)4 



(8) 
(9) 



in terms of parton matrix elements 

7^ 



^ . d'xY,m'M\Ph,s~i) 

-^Ph J s-l 



Tq{p,Ph) = Tr 

X {ph,S~l\i;,{x)\0)e'P-^ 



(10) 



Tg{l,Ph) 



d\Y.^Q\Al{x)\pu,S-l) 



S-l 



X (p,,5-l|A^(0)|0) 
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(11) 



They can be interpreted as single inclusive hadron multi- 
plicity distributions within fractional momentum Zh and 
Zh + dzh from parton fragmentation. Similarly, one can 
derive the momentum sum rules for dihadron fragmen- 
tation functions. Furthermore, one can also derive sum 



rules that relate dihadron fragmentation functions to sin- 
gle fragmentation functions. 

Given p_Lh ~ P±i + P±2 and q± = pj_i — pi_2, we have 
assumed that the total transverse momentum of the two 
hadrons p^_h is parallel to the parton's transverse momen- 
tum in the definition of parton dihadron fragmentation 
functions in Eqs. lj2J) and jSjl. Similarly, if we consider 
single hadron fragmentation, the hadron transverse mo- 
mentum p±^i is also parallel to the partons transverse 
momentum. We should then have the following identity. 



2^zIp+^ 



2-kzIp+'^ 



= 1, 



(12) 



which is just integration over the angle of the initial par- 
ton. Using the above identity, one can recast the mo- 
mentum integration as, 



dz2 ,2 d'^q±_ 
Z2 4(27r)'^ 



2dVi 



d^P2 



2E2{2^f 



(13) 



Here, the external quark momentum is p+ and hadrons 
momenta are Pi = zip^ and = Z2P^- One can also 
rewrite the (5-function S{zii — pf^/p'^) ~ (zi/ Zh)6{zi — 
Pi /P'^) the definition of the dihadron fragmentation 
functions. 

Since a dihadron fragmentation function is essentially 
a two-hadron multiplicity distribution, its sum rules must 
involve the information of the number of hadron pairs. 
An example of this is 

J2 J dzidz2D^^i;^-{zuZ2) = {N{N ~ 1)), (14) 

hi,h2 

where {N{N — 1)) is the second cumulant moment of the 
multiplicity distribution and 



{N) =Y, I dzDl^iz) 



(15) 



is the mean multiplicity from parton fragmentation. In- 
side the parton matrix elements of the fragmentation 
function, one can rewrite 



E 

h2,S-2 



E 



2-B2(2^)3 



\pi,P2,S -2){pi,p2,S -2\ 



d^P2 



2E2{2TTf 



ail ) (^^2 ) cih^ {P2 ) ah, {pi ) 



= 4,(piKi(pi)(^-i) 



(16) 
(17) 



where aj^ and dh are creation and annihilation operators 
for hadron h, and 

^ = E/^4(pK(p), (18) 
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is the total multiplicity operator. Here we have assumed 
hadrons as bosons. Using similar expression, 

J2\Pi.S'l){pi,S-l\=alipi)ahM, (19) 

S-l 

for single inclusive states, one can effectively have the 
identity, 

Using Eqs. l(T^ . ((T^ and 1^01) . one can obtain the follow- 
ing relationship, 

E / = (21) 

between dihadron and single hadron fragmentation func- 
tions. Using the momentum sum rules for the single 
hadron fragmentation functions, one has the following 
momentum sum rule for dihadron fragmentation func- 
tions, 

/ll./l2 ^ ' 

(22) 

Note that both (N) and (^(^-1)) in Eqs. ^ and ^ 
are not infrared safe and therefore not well defined in 
the collinear approximation for the fragmentation func- 
tions. Consequently, these sum rules for dihadron frag- 
mentation functions are also not well defined in collinear 
approximation. They, however, may provide useful phe- 
nomenological constraints for practical modeling of these 
fragmentation functions where collinear approximation 
might not be required. 

IV. DGLAP EVOLUTION 

One of the many successes of the factorized pQCD is 
the prediction of the scaling violation of the single inclu- 
sive fragmentation functions via the DGLAP equations. 
If the fragmentation functions are measured at an initial 
large scale Qo >> ^qcd, they can be predicted at any 
higher scale Q. In terms of the physical pictures afforded 
by the parton model, one imagines that the initiating 
quark loses its large virtuality through the radiation of 
multiple soft gluons and finally fragment into hadrons. 
Such a picture is validated within field theory by con- 
sidering the leading twist contributions of higher order 
diagrams in e+e^annihilation. The leading log (LL) con- 
tributions from all such diagrams are resummed and a set 
of coupled differential equations for the rate of change of 
the fragmentation functions with the momentum scale is 
derived. 



In Ref. 0, a QCD evolution equation for the non- 
singlet quark dihadron fragmentation function in an op- 
erator formalism was derived. The resulting evolution 
equation also admits a simple physical interpretation in 
terms of the diagrams shown in Figs. ^[5l|Hl The reader 
will note the new contribution from the diagram indi- 
cated in Fig. El the two detected hadrons emerging from 
the fragmentation of two separate partons. This contri- 
bution will turn out to be essential in determining the 
QCD evolution of dihadron fragmentation functions. 




FIG. 4: The contribution from the quark fragmentation to the 
NLO correction of the quark dihadron fragmentation func- 
tion. 




FIG. 5; The contribution of gluon fragmentation to the NLO 
correction of tire quark dihadron fragmentation function. 
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The resulting DGLAP evolution equation for the quark 
I 



91ogQ2 



Where D^{z,Q'^) and Dg{z,Q^) are the single in- 
clusive quark and gluon fragmentation functions, and 
Dg^''-^{zi,Z2,Q^) is the gluon dihadron fragmentation 
function [see Eq. ((SJ] . In the above equation Pq^qg (y) is 
the splitting function for a quark to radiate off a gluon 
and keep a fraction y of its initial forward light-cone mo- 
mentum. It is identical to the splitting function kernel of 
the DGLAP equation for the single fragmentation func- 
tion, i.e., 

P,^M = Cf{j:^^^. (24) 

The subscript indicates that the negative virtual cor- 
rection has been added within the splitting function. The 
splitting function in the second line of Eq. H23|) , Pq^qg , 
is identical to the above equation except that it lacks the 
negative virtual correction. The last splitting function is 
the probability for a quark to radiate off a gluon with a 
fraction y of its forward momentum and has an expres- 
sion identical to the case for the single fragmentation 
case. It should be pointed out in passing that Pq^gq{y) 
also has no virtual correction and thus no negative con- 
tribution. In the evolution of the dihadron fragmenta- 
tion function the sole negative contribution arises from 
the virtual piece in Pq^qg{y)] all other contributions are 
positive. Another difference from the case of the single 
inclusive fragmentation function are the measures l/y^ 
and l/y(l — y). These can be understood on the basis 
that the origin of the dihadron fragmentation function 
lies in the evaluation of a double differential inclusive 
cross section (see Sec. II of Ref. 0). We also point 
out that Eq. ^ may also be derived in the cut-vertex 
formalism via the renormalization of the bare cut- vertex 
presented in Fig. |21 

The evolution of the gluon dihadron fragmentation 
function may be derived in the cut-vertex formalism via 
the renormalization of the bare cut-vertex shown in Fig. 



dihadron fragmentation function is given as 



(23) 



I 

|31 It should be pointed out that the evolution equations 
for the quark and gluon dihadron fragmentation func- 
tions may be motivated using parton model arguments 




FIG. 6: The mixed contribution to the NLO correction of the 
quark dihadron fragmentation function. 



such as those of Refs. 0, l^^. The evolution equa- 
tion for the gluon dihadron fragmentation functions now 
includes four pieces: the gluon may split into a quark- 
antiquark pair or into two gluons. The detected hadrons 
may both emanate from the same parton (quark, anti- 
quark or gluon) or from the two different partons result- 
ing from the split, i.e., 



as 
27r 



+ 



1+Z2 y 

dy 
y(i-y) 

1-^^ dy 



^P (y)D'^^f^^ ( ^ ^ 



Pq^qg{y)D'^' 



Pq^qg{y)D'--{^,Q']D^g^ 



y 

Z2 



Z2 



i-y 



i-y 



dy 



^Pq^gq{y)D^^-'^- ( -,-,Q 



1+Z2 y 



Zl Z2 



y y 
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In the above equation, Uf is the number of flavors of 
quarks assumed. In the remainder of this paper we wiU 
assume n/ to be 3. Henceforth, the quark and antiquark 
fragmentation functions wiU not be distinguished; in both 
cases the singlet fragmentation function will be assumed 
i.e., Dq = Dq = {Dq + Dq)/2. The sphtting function 
Pg~tqq is the probability for the initial gluon to decay 
into a quark-antiquark pair with the quark carrying a 
fraction y of the forward momentum of the gluon. It is 
given as. 



g^gg 



(26) 



where Tp is the Casimir. It may be noted that inter- 
changing quark and antiquark may also be achieved by 
the switch y ^ 1 — y. In Eq. 1^3)) . Pg^gg{y) is the prob- 
ability for the gluon to split into two gluons, and has the 
expression 



Pg-.gg{y) = 2C A 



y 
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+ y(i - y) 



6 3 ^ 



(27) 



The presence of the '-|- '-function as well as the (5-function 
is the result of virtual contributions from gluon and 
quark-antiquark loops. The final splitting function in 
Eq. (|25|l is essentially the above splitting function with- 
out the virtual corrections, i.e.. 



P, 



g-*gg 



(y) - 2Ca 



y 



i-y 
y 



+ yi^-y) 



,(28) 



since there are no virtual corrections to the indepen- 
dent fragmentation contributions at the same order. As 
hadrons with finite momentum fractions originating from 
both gluons are detected, the range of values of the in- 
termediate momentum fraction y may approach neither 
nor 1. Therefore, there is no infrared divergence even 
without cancellation by virtual corrections. As one can 



see now there exist two new contributions to the evolu- 
tion of the gluon dihadron fragmentation function; from 
independent single fragmentation functions after a split 
to a quark-antiquark and to two gluons. Both these con- 
tributions are positive. This again will turn out to be an 
essential part in the equation that influences the QCD 
evolution of gluon dihadron fragmentation functions. 



V. RESULTS OF EVOLUTION: COMPARISON 
WITH EVENT GENERATORS 

To date, measurements of single inclusive cross sections 
in e+e^ collisions remain the primary set of data used 
in the parameterization of the fragmentation functions 
and testing their scaling violations. In many ways this 
provides an independent justification of the factorized 
pQCD approach to high energy collisions. The absence 
of any initial state interactions makes e~''e~experiments 
ideal for baseline measurements of parton fragmentation 
in vacuum, while allowing for tractable calculations of 
their evolution. However, no such measurements have 
been performed for double inclusive cross sections. In 
the absence of such measurements and in the interest of 
simplicity we turn to Monte Carlo event generators for 
both the extraction of the initial conditions and for com- 
parison to the numerical results of the evolution equation 
for dihadron fragmentation functions. 

Monte Carlo event generators such as JETSET 
have enjoyed great success as simulators of e+e~collision 
events. They provide reliable predictions not only for sin- 
gle inclusive measurements but also for many-particle ob- 
servables such as event shape and inter-jet particle flows. 
Hence, it is reasonable to assume that two particle cor- 
relations extracted from a "tuned" event generator will 
closely mimic such correlations measured in real experi- 
ments. The events generated will be restricted to two-jet 
events for a measurement of the quark dihadron frag- 
mentation function and to three-jet events for the gluon 
dihadron fragmentation function. Single inclusive hadron 
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fragmentation functions for both quarks and gluons will 
also be measured in the same set of events. 

In the extraction of the quark fragmentation functions 
three million dijet events distributed equally over three 
flavors were simulated using JETSET. The ^/s of the re- 
action was set to 20 GcV. In all such events the direction 
of the initiating quark and its initial virtuality were con- 
trolled. The virtuality of the initiating quark Q sets the 
scale of the fragmentation function. The forward light- 
cone momentum of the initiating quark was read off from 
the event list. In such experiments the single fragmenta- 
tion function is defined as 



1 dN{z) 



dz 



(29) 



In this paper, we restrict the flavor of the detected par- 
ticles to 7r+ and tt^. The variable z = p^/p'^, where 
Pq is the momentum of the fragmenting quark. The de- 
nominator Nevt is the number of events, whereas dN{z) 
represents the number of particles, over all events, that 
fall between z and z-\-dz of the momentum fraction. The 
dihadron fragmentation function is measured as 



Dq{zi, 22) = 



1 d2p(zi,Z2) 



dzidz2 



(30) 



where, d^ P{zi, Z2) represents the number of pairs of pi- 
ons with momentum fractions between (zi, zi +dzi) and 
{z2,Z2 + dz2). In constructing these functions, zi was 
always restricted to be larger than Z2- 




GLJON DIHADRON 



FIG. 7: The gluon dihadron fragmentation function. 

For the extraction of the gluon fragmentation function, 
an identical procedure as above was carried out with the 
sole restriction, that all events be three-jet events and a 
large fraction of the energy be concentrated in the gluon. 



The latter requirement guarantees that events are pre- 
dominantly composed of cases where the gluon and the 
quark antiquark pair are always contained in opposite 
hemispheres. Tracking the final hadrons in the gluon's 
hemisphere lead to the construction of the gluon frag- 
mentation function. It should be pointed out that in the 
Lund model of fragmentation, which is the underlying 
fragmentation model in JETSET, an out going gluon is 
represented by a kink in the fragmenting string which be- 
gins at the quark and terminates at the antiquark. Per- 
turbative showers within JETSET modify this picture 
and allow for multiple strings to form and fragment. A 
plot of the gluon dihadron fragmentation function thus 
extracted is shown in Fig. [7| 




FIG. 8: The result of At of the function of Eq. 1^3 to 
Fig.Q 

For the convenience of providing initial conditions to 
the evolution equations, the fragmentation function of 
Fig. [7| is parametrized by fitting to a function of the 
type: 



D{zi, Z2) 



Nz^'z^^izi 



\P2 



(1-zi 



Z2r'{l~Z^f'{l-Z2) 
Z2f'. (31) 



There are seven parameters in the above fit: 
iV, Q!i, a2, as, /?!, /?3. The reader will note that 
the structure of the function is a simple generalization of 
the parameterization of the single fragmentation func- 
tions. The fit is carried out by the method of minimum 
logarithm of . The logarithm ensures that the fit is 
better at larger values of zi and 22. A plot of the final 
result of the fit for the gluon dihadron fragmentation 
function is shown in Fig. |S1 One notes that the fit 
function of Eq. (|31|l mimics the function closely except 
at very small zi and 22- This is unimportant as in the 
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evolution we will always restrict our attention away 
from very small Zi and Z2. We will not explicitly show 
the quark dihadron fragmentation function and its fit 
here. Suffice to say that the fit is even better for the 
quark, which, as is the case for the single fragmentation 
function, generically has a harder spectrum in momen- 
tum fraction. Values for the various parameters of the 
fits to the quark and gluon fragmentation functions at 
the two different values of the used in this paper are 
presented in tabled 

One unique feature of the DGLAP evolution equa- 
tions for dihadron fragmentation functions is that the 
equations couple dihadron to single hadron fragmenta- 
tion functions. The single hadron fragmentation func- 
tions themselves evolve independently according to their 
own DGLAP evolution equations. We find that the single 
hadron fragmentation functions obtained from JETSET 
simulations can be described very well by Binneswies- 
Kniehl-Kramer(BKK) parameterization [23 of the actual 
experimental data. Therefore, in our numerical study of 
the evolution equations for dihadron fragmentation func- 
tions, we will simply use the BKK parameterization for 
the single hadron fragmentation functions and their evo- 
lution with the momentum scale Q. 

The fit function shown in Fig. |H1 as well as its ana- 
logue for the quark will provide the initial conditions to 
the differential equations outlined in Eqs. (|25|l and 123|l 
respectively. In both cases, the initial scale of the frag- 
mentation functions is set to Qq = 2 GeV^. This cor- 
responds to log((5o) = 0.693. These are shown as the 
filled triangles for the quarks and the filled squares for 
the gluons in Fig|5| We have chosen a fixed zi = 0.5 and 
let Z2 vary from to 1 — zi, since several experiments 
that measure medium modification of the dihadron frag- 
mentation functions, which we will discuss in separate 
studies, have similar kinematic range. Note the orders 
of magnitude difference between the two distributions. 
This feature holds for most of the range of zi and Z2 ex- 
cept at very small momentum fractions where the gluon 
fragmentation function overtakes that of the quark. 

Similar to the procedure carried out in Ref. Q for 
non-singlet quark dihadron fragmentation functions, re- 
sults of the evolution will be presented in increments of 
log((5^) = 1- As expected, we note a softening of the 
spectrum with rising scale. We terminate the evolution 
at log{Q^) = 4.693, corresponding to scale = 109 
GeV^. To compare the results of the evolution equa- 
tions in Eas. H23|) and H25|l to Monte Carlo event simu- 
lations, dihadron fragmentation functions at the highest 
scale are once again extracted from JETSET. The num- 
ber of events used to sample the fragmentation functions 
and the method of the fit remain identical to the case at 
lower described above. The results are presented as 
open triangles for the quarks and open squares for the 
gluons. As shown they are in excellent agreement with 
the results of the evolution equations. This provides the 
most crucial test of the evolution equations we have pre- 
sented in this paper. 
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FIG. 9: Results of the evolution of the quark and gluon 
fragmentation function {Dq{zi, Z2), Dg{zi, Z2)). In all cases 
21 is held fixed at 0.5. 



Parton 


TV 


Ql 


02 


as 


Pi 


/32 


/33 


quark 
2 GeV^ 


4.080 


-0.673 


-0.440 


-0.707 


0.196 


1.717 


1.359 


quark 
109 GeV^ 


5.872 


-1.103 


-0.425 


-0.436 


0.410 


2.997 


2.164 


gluon 
2 GeV^ 


8.000 


-6.246 


-1.319 


6.736 


4.324 


15.214 


1.351 


gluon 
109 GeV^ 


1.090 


-8.848 


-1.430 


8.568 


6.216 


22.031 


-0.145 



TABLE I: Values of different parameters used in the fit to 
the fragmentation functions. 



In the recent experimental studies e.g., two high px 
particle production in p+p, p + A and A + A collisions at 
RHIC, or two particle correlation in inclusive DIS by the 
HERMES experiment, one usually measures dihadron 
correlations in the form of inclusive spectra of associated 
particles produced in correlation with a high momentum 
trigger hadron. These measurements are essentially the 
ratios of the number of trigger and associated particle 
pairs divided by the number of triggers. They are equiva- 
lent to the ratios of the dihadron fragmentation functions 
to the single fragmentation functions. This quantity is 
plotted in Fig. ^|for the same range of scales as that in 
FigEl The value of zi is once again held fixed at 0.5. The 
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FIG. 10: Results of the evolution of the ratio of the dihadron 
fragmentation function to the single fragmentation function 
for quarks and gluons. 



single hadron fragmentation functions are given by BKK 
parameterization at same which can also be de- 
rived from the same Monte Carlo simulations or evolved 
with the DGLAP equations for single hadron fragmenta- 
tion functions. We have checked that all three methods 
give almost identical results. While it may come as no 
surprise that the evolution of this ratio is also predicted 
by the dihadron DGLAP evolution equations, it must 
be pointed out that the ratio of the dihadron and single 
hadron fragmentation functions do not display substan- 
tial change with variation of the scale. This is especially 
true for the gluon fragmentation function. 



VI. DISCUSSIONS AND CONCLUSIONS 

In this study, a very general discussion of the properties 
of the dihadron fragmentation functions has been carried 
out. Various results have been outlined: the momentum 
sum rules have been discussed, evolution equations for 
both the singlet quark and gluon dihadron fragmenta- 
tion functions have been derived and solved numerically, 
comparisons with results from Monte Carlo event gen- 
erators (in the absence of experimental measurements) 
are made. Very good agreement has been obtained from 
such comparisons. In the remaining, we cast a backward 
glance at dihadron fragmentation functions, summariz- 



ing the sentinel points of the aforementioned analysis, 
and present an overview of future studies. 

Factorization ensures that fragmentation functions are 
well defined and universal. Once defined and measured in 
a given process they must admit the same definition and 
hence numerical value in another experiment. Such def- 
initions in light-cone gauge were presented in Ref. in 
LO of QCD. Gauge invariance of the definition is appar- 
ent and the generalization involves little more than the 
similar procedure invoked in the case of the single inclu- 
sive fragmentation function P|: The partonic operators 
at two space time points were connected by a Wilson line 
of the gauge field. In the case of the gluon fragmentation 
function, the gluon vector potential was replaced with 
gauge invariant field operators. 

Single inclusive fragmentation functions, which have 
the interpretation of single inclusive hadron distributions, 
obey well defined sum rules. The momentum sum rules 
for the dihadron fragmentation functions, which have the 
interpretation of the pair multiplicity in a jet, have been 
derived in Sec. III. One can relate the dihadron to sin- 
gle hadron fragmentation functions through these sum 
rules. However, they will involve first and second order 
cumulant moments of the multiplicity distribution since 
dihadron fragmentation functions involve pair multiplic- 
ity. The appearance of the cumulant moments makes 
the sum rules for dihadron fragmentation functions less 
rigorous since they are not well defined in the coUincar 
factorization approximation. However, in practice, they 
will provide useful guidance for phenomenological mod- 
eling of the dihadron fragmentation functions. 

Though fragmentation functions are well defined in 
terms of parton matrix elements, yet they still involve 
non-perturbative objects. Therefore, the exact form of 
a fragmentation function may not be estimated purely 
from QCD. They, however, can be measured at a given 
energy scale for a whole set of continuous forward mo- 
mentum fractions. Once such a measurement is made, 
one may derive the value at any other higher scale via 
the evolution equations in pQCD. Comparing these with 
experimental results at higher scales amounts to a test 
of the QCD evolution equations. In this paper, the lack 
of experimental results was circumvented by the use of 
the JETSET event generator which is tuned to fit a vari- 
ety of experimental data on e+e~annihilation. Note that 
event generators can reproduce not only the single inclu- 
sive measurements in e"'"e~events, but also a multitude of 
many particle observables. We used dihadron fragmen- 
tation functions from JETSET simulations at the scale 
of Qo = 2 GeV^ as the initial condition for the DGLAP 
evolution equations. We demonstrated that the results 
at higher scales, e.g., = 109 GeV^ from the numerical 
solution to the evolution equations are in excellent agree- 
ment with that extracted from JETSET simulations at 
the same scale, as shown in Figs. O and ^| This com- 
parison provided a stringent test to the new evolution 
equations presented in this study. 

The effect and validity of each component of the 



DGLAP equations may be further tested by a simple 
exercise. The various terms presented in Eqs. (|23|l and 
(j25(l may be broadly divided into two categories. The 
regular components, which we call correlated two hadron 
fragmentation, are contained in the first and last lines of 
Eq. (123) and the first and fourth lines of Eq. . These 
are simple generalizations of the single inclusive DGLAP 
equations and depend solely on dihadron fragmentation 
functions. The new components, which we call indepen- 
dent two hadron fragmentation, make up the second and 
third line of Eq. (|23|) and the second, third and last line 
of Eq. H25|) . These depend on products of single inclu- 
sive fragmentation functions, and therefore couple the 
evolution of the dihadron fragmentation functions with 
that of the single fragmentation functions. We may com- 
pute the evolution of dihadron fragmentation functions 
without the contribution from the components of inde- 
pendent fragmentation. This leads to the results pre- 
sented in Fig. As the reader may note, the calcu- 
lated evolution is now a poorer fit to the results from 
JETSET. While this may not be very noticeable for the 
quark fragmentation function, the gluon fragmentation 
function which is at least an order of magnitude smaller 
than the quark fragmentation function shows a very no- 
ticeable difference. The evolution of the gluon fragmen- 
tation function depends on two separate pieces of inde- 
pendent fragmentation. One receives contributions from 
the product of two gluon single fragmentation functions, 
while another from the product of a quark and an an- 
tiquark distribution functions. It is the removal of lat- 
ter that produces the majority of the difference in the 
evolution of the gluon dihadron fragmentation function 
between Fig. and Fig. ^] This is due to the fact that 
we restrict our attention to large z where the quark frag- 
mentation function dominates over the gluon, and the 
product of two single quark fragmentation functions is 
larger than the gluon dihadron fragmentation function 
at large zi and Z2. Thus each component of the evolu- 
tion has a role to play in the change of the fragmentation 
function from one scale to another. The contributions 
from the new components of independent fragmentation 
are always positive. Therefore, they slow down the scale 
evolution of the dihadron fragmentation functions, par- 
ticularly for gluon jets at large zi and Z2. This is also 
the reason why the triggered distributions (the ratio be- 
tween dihadron and single hadron fragmentation func- 
tions) change very little with the scale for gluon jets as 
shown in Fig. |51 

With the question of evolution of the dihadron frag- 
mentation functions set aside, we now focus on certain 
general properties of the fragmentation function as has 
been exposed by this analysis. It may be noted from 
Fig. ^1 that the ratio of the dihadron fragmentation 
function to the single fragmentation function of the lead- 
ing hadron (D^i''^ (zi, za, Q^)/!)^! (zi, Q^) ) shows little 
change as a function of even as — s- 100 GeV^. This 
is especially true of the ratio in the case of the gluon frag- 
mentation function, which shows practically no change 
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FIG. 11: Results of the DGLAP evolution of the dihadron 
fragmentation function without the independent fragmenta- 
tion terms. AH parameters are the same as Fig. |^ 



with over the range of zi and Z2 explored. The re- 
sults of evolution are strongly dependent, however, on 
the initial conditions and thus on the actual values of zi 
and Z2- This is consistent with the observations noted 
in our previous study of quark non-singlet dihadron frag- 
mentation functions 0- 

In recent experiments at RHIC 3, correlations of two 
high p± hadrons have been measured in p+p, d + Au and 
Au + Au collisions. Hadrons with p± > 4 GeV are used 
as a trigger. Once such a "leading" particle is identi- 
fied, the experiment measures the differential probability 
of the "next-to-leading" hadrons or associated high p± 
particles emanating from the same collision with trans- 
verse momentum in the range 2 < p± < 4 GeV at a 
given azimuthal angle (j) with respect to the direction of 
the leading particle. In Ref. 0, results for dN/d(j)/Ntrig 
are measured. While a large suppression was noted for 
w TT (away-side) in central collisions, almost no change 
was seen in the vicinity of ~ (near side). Assuming 
both hadrons (leading and associated) with the same- 
side correlation come from fragmentation of a jet, the 
initial jet energy is at least about 7 to 10 GeV. This is 
about the same range of momentum scale we have consid- 
ered in studying the evolution of dihadron fragmentation 
functions. This range of is also relevant to the exper- 
iments of DIS on nuclei Q , where medium modification 
of the ratio of the dihadron fragmentation function to the 
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single fragmentation function of the leading hadron was 
reported. A double ratio of the above quantity Ra/ Rd 
with a large nucleus versus that in deuteron showed min- 
imal change with the atomic number. A complete un- 
derstanding of the relevance of this observation requires 
a repeat of the calculation presented in this article with 
the inclusion of a medium modification. While the study 
with medium modification will be presented in a forth- 
coming article, we can already see the trend by analyzing 
the analogous effects of QCD evolution in the vacuum. 

If the results of Ref. Q for dN/d(j)/Ntrig are integrated 
over the angle in the vicinity of = 0, one should 
obtain, in effect, the ratio of the dihadron fragmentation 
function to the single fragmentation function. This is, 
no doubt, based on the assumption that both the leading 
and the next-to-leading particles emerge from the same 
parent parton. In the results of Ref. Q the ratio is 
essentially integrated from Z2 ~ 0.3, zi sa 0.6 to Z2 ~ 
0.4, zi « 0.4. Within this kinematic range, there is not 
much change of the ratio of the fragmentation functions 
as a function of the due to gluon bremsstrahlung in 
vacuum. Therefore, it may not be entirely surprising 
that no variation was noted in the same side two hadron 
correlations in Ref. due to medium induced gluon 
bremsstrahlung, in particular if one takes into account 
the trigger bias caused by parton energy loss [2^ . 

It has been pointed out recently that the ratio 
of the dihadron fragmentation function to the sin- 
gle fragmentation function of the leading hadron (i.e., 
D{zi, Z2)/ D{zi)) may be numerically similar to the sin- 
gle fra gme ntation function of the associated hadron i.e., 
D(z2) |25j. This was noted experimentally and in simu- 
lations. We point out that there is indeed some truth 
to this observation. In Fig. E| we plot the ratios of 
the dihadron fragmentation function to the fragmenta- 
tion function of the leading hadron for the quark and 
gluon. This is compared with the single inclusive frag- 
mentation functions of a quark (solid lines) and gluon 
(dashed lines). In the range < Z2 < 0.3, the single 
inclusive fragmentation functions do indeed closely ap- 
proximate the ratio Dq{zi, Z2) / Dq{zi). This fact is how- 
ever only true for the quark; no such similarity is noted 
for the gluon. The primary reason for the difference be- 
tween Dq{z2) and Dq{zi, Z2) / Dq{zi) for Z2 > 0.3 is the 
differences in kinematic bounds experienced by the two 
quantities: Dq{z2) ^ as Z2 — > 1 whereas Dq{zi, Z2) — > 
as 2:2 ^ 1 — zi. In Fig. ll2[ zi is held fixed at 0.5, thus the 
maximum value of Z2 = 0.5. This differences in kinematic 
bounds may be circumvented by considering a rescaled 
single fragmentation function (22/(1 — 2i))/(l — zi) = 
2Dq(2z2) which experiences a similar kinematic bound as 
Dq{zi, Z2) / Dq{zi) for zi = 0.5. This rescaled function is 
represented by the dot-dashed lines in Fig. E| which is 
very similar to Dq{zi, Z2) / Dq{zi) with the same momen- 
tum scale Q. 

Finally, we point out that in the derivation of results 
presented in this paper, we depended on the assumption 
that the energy scales of the processes in question were 
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FIG. 12: Comparisons between the ratio of the dihadron 
fragmentation function to the single fragmentation functions 
of the leading hadron (indicated as D/S) and the fragmenta- 
tion function of the associated hadron (indicated as S). The 
D/S curves are the same as in Fig. 1101 The single fragmenta- 
tion functions are obtained from the BKK parameterization of 
Ref. • The dashed lines are for the gluon and the solid lines 
for the quark fragmentation function. The dot-dashed lines 
represent a rescaled quark fragmentation function 2Dg^ (2^2). 
The of a particular dot-dashed line is the same as the 
of the solid line in the same order. 



high enough for the applicability of pQCD methods. We 
also required the fragmentation functions to be defined 
at a scale larger than a semihard scale n±, such that 
^QCD << A^i Q^- This semi-hard scale fi± restricts 
the relative transverse momentum between hadrons such 
that they are considered as from one parton jet. We thus 
required to be high enough for a hierarchy of scales 
to exist. The evolution equations are applicable solely 
when this hierarchy of scales exists. 
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